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ABSTRACT 

We report the discovery of planets orbiting two bright, nearby early K dwarf stars, HD 97658 and 
Gl 785. These planets were detected by Keplerian modelling of radial velocities measured with Keck- 
HIRES for the NASA-UC Eta-Earth Survey. HD 97658 b is a close-in super-Earth with minimum 
mass Msini = 8.2 ± 1.2 M®, orbital period P = 9.494 ± 0.005 d, and an orbit that is consistent 
with circular. G1785b is a Neptune-mass planet with Msini = 21.6 ± 2.0 M e , P = 74.39 ± 0.12d, 
and orbital eccentricity e = 0.30 ± 0.09. Photometric observations with the T12 0.8 m automatic 
photometric telescope at Fairborn Observatory show that HD 97658 is photometrically constant at 
the radial velocity period to 0.09 mmag, supporting the existence of the planet. 

Subject headings: planetary systems — stars: individual (HD 97658, G1785) — techniques: radial 
velocity 



1. INTRODUCTION 

Radial velocity (RV) searches for extrasolar planets 
are discovering less massive planets by taking advan- 
tage of improved instrumental precision, higher observa- 
tional cadence, and diagnostics to identify spurious sig- 
nals. These discoveries include pla nets with minimum 
masses (Msini) as low as 1.9 Mq ([Mayor et al.l l2009f) 
and sy stems of mult i ple lo w -mass plane t s (lLovis et al.l 
[20061: iFischer et all 120081 IVoet et al.1 12010^ To 
date, 15 planets with Msini<lOM0 and 18 plan- 
ets with Msini = 10-30 M m have be en discovered by 
the RV technique (jWright et al.l 120101 Exoplanet Orbit 
DatabaseFT) . Transiting searches for extrasolar planets 
have detected Neptun e-mass planets dBakos e t al.ll2010t 
Hartman et al.l l2010fl and sup er-Earths (jLeger et al.1 
2009t iCharbonneau et al.l 120091) . The initial data re- 
lease from the Kepler mission shows substantially in- 
creasing planet occur rence with decreasing planet radius 
(jBorucki et al.l 120101 ) . Using the large number of low- 
mass planets, we can statistically study planet proper- 
ties, occurrence rates, and parameter correlations in ways 
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23l^% for planets in the mass range 0.5- 



previously only possible with higher mass gas-giant plan- 
ets. 

We recently completed an analysis of close-in planet 
occurrence for 166 G- and K-type dw arf stars in the Eta- 
Earth Survey (jHoward et al.l l2010d) . We studied the 
planet detections and non-detections on a star-by-star 
basis, estimating search completeness. We detected in- 
creasing planet occurrence with decreasing planet mass 
over the mass range 3-1000 M ffi for planets with orbital 
periods P < 50 d. We parameterized the planet mass dis- 
tribution with a power law model from which we extrap- 
olated the occurrence rate of close-in Earth-mass planets, 
giving ry e = 23±igS 
2.0M ffi with P < 50 d. 

Our study also addressed a key predictio n of popu- 
latio n synthesis models of p lanet formation ()Ida fc Linl 
12004 120081 : iMordasini et al.ll2009| )— the expected dearth 
of close-in, low-mass planets. The "desert" emerges 
in the simulations from fast migration and accelerating 
planet growth. Most planets are born near or beyond the 
ice line and those that grow to a critical mass of several 
Earth masses either rapidly spiral inward to the host star 
or undergo runaway gas accretion and become massive 
gas giants. Our measurements contradict this prediction; 
we found the highest occurrence rate for planets where 
theory predicted a dearth, in the regime of 3-30 M e and 
P < 50 d. Population synthesis models of planet forma- 
tion are currently unable to explain the distribution of 
low-mass planets. 

To measure the planet occurrence rate as a function 
of planet mass, our study included previously detected 
planets, as well as unannounced "planet candidates" 
(|Howard~et al. 2010c). Including candidates was neces- 
sary to reliably estimate occurrence fractions for low- 
mass planets, even though the candidates had formal 
false alarm probabilities (FAPs) as large as 5% at the 
time of our analysis (June 2010). Such an FAP implies 
that the planet is very likely to exist, but it's too high for 
the secure announcement of a definite planet detection 
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with well-measured orbital parameters. Since then, we 
continued to intensively observe the planet candidates. 
Based on the new confirmatory data we report two of 
them here as bona fide planets. We present HD 97658 b, 
a close-in, super-Earth p lanet identified as "Candidate 3" 
in Howar d et al.1 (|2010cf ). and G1785b, a Neptune-mass 
planet identified as "Candidate 7". 

Below we describe the host stars (Section [2]) and the 
RV measurements (Section [3]) . We analyze these mea- 
surements with Keplerian models and assess the prob- 
ability of spurious detections by computing false alarm 
probabilities (Sections @] and [S]). We describe photomet- 
ric observations of HD 97658 and the limits they impose 
on planetary transits (Section [6]). We discuss the radii 
of these planets and a trend in the host star metallicities 
among low-mass planets (Section [7]) . 

2. STELLAR PROPERTIES 

We used Spectroscopy Made Easy (IValenti fe Piskunovl 
119960 to fit high-resolution spectra of HD 97658 
(HIP 54906, GJ3651) and G1785 (HD 192310, 
HIP 99825), using the wavelength inte r vals, line 
data, and methodology of Valenti fc Fischer (200a). We 
further constrained surface gra yity using Yonsei-Yale 
(Y 2 ) stellar structure models (De marque et al.l 120041) 
and revised Hipparcos parall axes (Ivan Leeuwenl I2007D . 
using the iterative method of IValenti et al.l (|2009f ) . The 
resulting stellar parameters listed in Table Q] are effective 
temperature, surface gravity, iron abundance, projected 
rotational velocity, mass, radius, and luminosity. Both 
stars are K dwarfs on the main sequence. 

HD 97658 lies 0.46 mag below the Hipparcos average 
main sequence (My versus B — V) as defined by [Wright 
(2005) . This location is consistent with the low metallic- 
ity of [Fe/H] = -0.23 ± 0.03. G1785 is 0.06 mag above 
the Hipparcos average main sequence, consistent with its 
slightly super-solar metallicity of [Fe/H] = +0.08 ± 0.03. 

Measurements of the cores of the Ca II H & K lines of 
each spectrum show low levels of chromospheric activity, 
as quantified by the Shk and log i? H K . These chromo- 
spheric indices show long-term trends over the six years 
of measurements, possibly partial activity cycles, so we 
list ranges of activity indices in Table [TJ We detect a 
weak correlation between individual RVs and Shk mea- 
surements for HD 97658, but not for G1785. This cor- 
relation, with a Pearson linear correlation coefficient of 
r = +0.35, docs not appear to affect the Keplerian fit 
of HD 97658 b because the Shk time series has negligible 
Fourier power at or near the adopted orbital period, even 
when the lo ng-term activity tren d is re moved. 

Following llsaacson fe Fischerl (|2010l ). and based on 
Shk, My, and B — V, we estimate an RV jitter of 
1.5ms" 1 for these stars. This empirical estimate is based 
on an ensemble of stars with similar characteristics and 
accounts for RV variability due to rotational modulation 
of stellar surface features, stellar pulsation, undetected 
planets, and u ncorrected systematic errors in the veloc- 
ity reduction (|Saar et alJll998t lWrightll2005] ). Jitter is 
added in quadrature to the RV measurement uncertain- 
ties for Keplerian modelling. 

3. KECK-HIRES VELOCITY MEASUREMENTS 

We observed HD 9 7658 and G1785 with the HIRES 
echelle spectrometer (|Vogt et al.lll994l ) on the 10-m Keck 
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Fig. 1. — Keck-HIRES spectra of the Ca II H line of the early 
K dwarf stars HD 97658 and G1785. Slight line core emission near 
3968 A indicates modest chromospheric activity. 

I telescope. The observations of each star span six years 
(2004-2010). All observations were made with an io- 
dine cell mounted directly in front of the spectrometer 
entrance slit. The dense set of molecular absorption 
lines imprinted on the stellar spectra provide a robust 
wavelength fiducial against which Doppler shifts are mea- 
sured, as well as strong constraints on the shape of the 
spectrome ter instrumental profil e at the time of each ob- 
servation (|Marcv fe Butledll99HValenti et al.lfl995l) . 

We measured the Doppler shift of each star-times- 
iodine spectrum u s ing a modelling proced ure descended 
from IButler et al.l (|1996|) as described in IHoward et al.l 
(|2010dD . The velocity and corresponding uncertainty for 
each observation is based on separate measurements for 
~700 spectral chunks each 2 A wide. Once the two plan- 
ets announced here emerged as candidates (about two 
years ago) we increased the nightly cadence of measure- 
ments and made three consecutive observations per night 
to reduce the Poisson noise from photon statistics. We 
calculate one mean velocity for multiple observations in 
a 2 hr interval. 

The highest RV measurement precision using Keck- 
HIRES has been achieved on chromospherically inactive 
late G and early K dwarfs, like the two stars presented 
here. The quietest of these stars are stable over many 
years at the ~l.5-2.0ms" 1 level (jHoward et al.l 120091 
2010a jJ; velocity residuals are due to astrophysical per- 
turbations, instrumental/systematic errors, and Poisson 
noise. All of the measurements reported here were made 
after the HIRES CCD upgrade in 2004 August and do 
not suffer from the higher noise and systematic errors 
that limited the precision of pre-upgrade measurements 
to ~2-3ms _1 for most stars. 

For each star we constructed a single-planet Keple- 
rian model using the o rbit fitting techniques described in 
IHoward et al.l (|2010af) and the partial ly linearized, least- 
square s fitting procedure described in (Wright fc Howard! 
(2009). The Keplerian parameter uncertainties for 
each planet were d erived using a Monte Carlo method 
(|Marcv et al.l I200"5l ) and do not account for correlations 
between parameter errors. Uncertainties in Msini re- 
flect uncertainties in M* and the orbital parameters. 

4. HD 97658 
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TABLE 1 

Stellar. Properties of HD 97658 and Gl 785 



Parameter 


HD 97658 


G1785 


r> *f v 1 fir c\ 
O L L 1 ctl LVjJU 


Kl V 


Kl V 


Mir 


6.27 


6 13 


B — V 


0.80 


0.78 


v 


7.78 


5.73 


J 


6.203 


4.112 


H 


5.821 


3.582 


K 


5.734 


3.501 


Distance (pc) 


21.1 ±0.33 


8.911 ± 0.024 


T cff (K) 


5170 ± 44 


5144 ± 50 


log 9 


4.63 ±0.06 


4.60 ±0.06 


[Fe/H] 


-0.23 ± 0.03 


±0.08 ± 0.03 


vs'mi (kms ) 


0.5 ± 0.5 


0.5 ±0.5 


L* (L ) 


0.34 ± 0.02 


0.30 ±0.02 


AT* (M ) 


0.85 ± 0.02 


0.78 ± 0.02 


R* (Re) 


0.73 ± 0.02 


0.68 ± 0.02 




-4.95 to -5.00 


-4.90 to -5.02 


Shk 


0.169 to 0.197 


0.169 to 0.226 




100 
Period (d) 

Fig. 2. — Lomb-Scargle periodogram of RV measurements of 
HD 97658. The tall peak near P = 9.494 d suggests a planet with 
that orbital period. 



TABLE 2 
Orbital Solution for HD 97658 b 

Parameter Value 



10 



P (days) 9.494 ± 0.005 

T c (JD - 2,440,000) 15375.01 ± 0.64 

e<1 =°-° o 

K (ms" 1 ) 2.75 ±0.39 .O 

Msini (M e ) 8.2 ± 1.2 3 

a (AU) 0.0831 ±0.0011 

N a hs (binned) 96 

Median binned uncertainty (ms — 1 ) 0.74 

Assumed jitter (ms~ 1 ) 1.5 —10 

a (ms -1 ) 2.78 

v/xf 1-59 



a We adopt a circular orbital solution for this planet. 

The RVs and Shk values from Keck-HIRES are listed 
in Table |H Figure H| show s a Lomb-Scargle periodogram 
(|Lomblll97 6; Scarglc 1982}) of the RVs with a substantial 
peak at 9.494 d. We used that period, as well as a wide 
variety of other tr ial periods, as seeds for the Keplerian 
fitting algorithm (j Wright fc Howard I2009D . Our search 
identified the single-planet orbital solution listed in Table 
[2] as the best fit. 

We also tried fitting the RVs with an eccentric Keple- 
rian model and found a best-fit solution with a nearly 
identical orbital period and e = 0.17 ± 0.17, which 
is consistent with circular at the 1-a level. The de- 
tection of nonzero eccentricity with better than 95% 
confidence (2-c) requires approximately e/a e > 2.45, 
where a e = a/K ■ (2/N) - 5 , a is the measurement un- 
certainty (including jitter), and N is t he number of 
uniformly phase-distributed observations (|Valenti et al.l 
l2009tlLucv fc Sweenevlll97l"l) . Our measurements do not 
meet this criterion. Furthermore, the eccentric model 
does not improve \/xJ from the circular model. We 
adopt the circular orbit model in Table [5J 

We considered the null hypothesis — that the observed 
RVs are the chance arrangement of random velocities 
masquerading as a coherent signal — by calculating two 
false alarm probabilities (FAPs). Using the method de- 
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Orbital Phase 

Fig. 3. — Single-planet model for the RVs of HD 97658, as mea- 
sured by Keck-HIRES. The dashed line shows the best-fit circu- 
lar orbital solution. Filled circles represent phased measurements 
while the open circles represent the same velocities wrapped one 
orbital phase. The error bars show the quadrature sum of mea- 
surement uncertainties and 1.5ms" 1 jitter. 

scribed in iHoward et al.l ()2010al ) , we computed the im- 
provement in A\ 2 from a constant velocity model to a 
Keplerian model for 10 3 scrambled data sets. In the first 
FAP test, we allowed for eccentric single-planet orbital 
solutions in the scrambled data sets. We found that three 
scrambled data sets had a larger Ax 2 than the measured 
velocities, implying an FAP of ^0.003 for this scenario. 
When we restricted the search for orbital solutions to cir- 
cular orbits, none of the scrambled data sets had a larger 
A\ 2 than measured velocities, implying an FAP of less 
than -0.001. 

The rms of 2.78 ms -1 about the single-planet model 
is relatively high compared to our adopted jitter of 
1.5 ms -1 for this chromospherically quiet K dwarf star. 
This suggests that the measured RVs are compatible 
with additional detectable planets. We computed a pe- 
riodogram of the RV residuals to the single-planet fit 
and found several periods with considerable power in the 
range ~40-200d. These peaks correspond to Doppler 
signals with ~1.5-3ms _1 semiamplitudes. We consid- 
ered two-planet orbital solutions with Pb seeded with the 
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TABLE 3 
Orbital Solution for. Gl785b 



Parameter 



Value 



P (days) 

T c (JD - 2,440,000) 

e 

K (ins" 1 ) 
Msini (M®) 
a (AU) 

N bs (binned) 

Median binned uncertainty (ms — 1 ) 
Assumed jitter (ms -1 ) 
<x (ms -1 ) 



74.39 ±0.12 
15173.2 ± 2.0 
0.30 ± 0.09 
4.07 ±0.41 
21.6 ±2.0 
0.319 ± 0.005 
73 
0.68 
1.5 
2.06 
1.17 



1 




100 
Period (d) 

Fig. 4. — Lomb-Scargle periodogram of RV measurements of 
G1785. The tall peak near P = 74.4 d suggests a planet with 
that orbital period. 

best-fit value from the single-planet model and P c seeded 
with peaks in the residual periodogram. We allowed all 
orbital parameters includin g eccentricities to float i n the 
two-planet fitting process (jWright fc How ard 2009[). No 
two-planet solutions were found with an FAP below 5%. 
We will continue to observe this star in search of addi- 
tional planets. 

5. GL785 

The RVs and Shk values from the Keck-HIRES mea- 
surements of Gl 785 are listed in Table [5] Figure [4] shows 
a Lomb-Scargle periodogram (lLonibl [T976: Scargle 1982) 
of the RVs with a substantial peak near 74.4 d. We 
identify the peaks near 1.0 d as stroboscopic aliases of 
the sider eal day with the 74.4 d sig nal and other long 
periods ([Dawson & Fabrvckvl 120101) . We used 74.4 d, 
as well as a wide variety of other periods, as seed pe- 
riods for the single-plan et Keplerian fitting algorithm 
(jWright fc Howard! 12009ft . Our search identified the 
single-planet orbital solution listed in Table [2] as the best 
fit. The orbital eccentricity of 0.30 ± 0.09 is significant 
at the 3-<7 level. 

We considered the null hypothesis for the observed pe- 
riodic signal in the measured RVs of Gl 785 by computing 
an FAP using the method described in Section |U includ- 
ing allowing for eccentric solutions with the scrambled 
data sets. We found that none of the 10 3 scrambled data 
sets had a larger A% 2 than the measured velocities, im- 
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Fig. 5. — Single-planet model for the RVs of G1785, as measured 
by Keck-HIRES. The dashed line shows the best-fit eccentric or- 
bital solution. Symbols have the same meanings as in Figure [3] 

plying an FAP of less than ~0.001. 

With an rms of 2.06 ms" 1 and a featureless peri- 
odogram of velocity residuals to the one-planet model, 
we do not see evidence for a second detectable planet 
orbiting Gl 785. 

6. PHOTOMETRIC OBSERVATIONS 

We also acquired photometric observations of 
HD 97658 with the T12 0.80 m automatic photometric 
telescope (APT) , one of several automatic telescopes op- 
erated by Te nnessee State Univ ersity (TSU) at Fairborn 
Observatory ([Eaton et al.ll2003D . G1785 is too far South 
for this observatory. The APTs can detect short-term, 
low-amplitude brightness changes in solar-type stars 
resulting from rotational modulation in the visibility 
of active regions, such as starspots and plages (e.g., 
iHenrv et all ll995b[ ) and can also detect longer-term 
variations produced by the growth and decay of in- 
dividual active regions and the occurrenc e of stellar 
magn etic cycles (e.g., IHenrv et al.l Il995at IHall et ahl 
[20091) . The TSU APTs can disprove the hypothesis 
that RV variations are caused by s tellar activity, rathe r 
than planetary reflex motion (e.g.. IHenrv et al.l [2000a) . 
Several cases of apparent periodic RV variations in 
solar-type stars induced by the pre sence of photo s pheric 
stars pots have bee n discu ssed by IQueloz et al.l ([2001D 
and IPaulson et al.l (|2004[) . Photometry of planetary 
candidate host stars is also useful t o search for transits 
of the planetary companions (e.g., IHenrv et al.l l2000bl : 
ISato et al.ll2005Tl . 

The T12 0.80 m APT is equipped with a two-channel 
photometer that uses two EMI 9124QB bi-alkali photo- 
multiplier tubes (PMTs) to make simultaneous measure- 
ments of a star in the Stromgren b and y passbands. 
The T12 A PT is function ally identical to the T8 APT 
described in IHenrv! (|1999l ). The final data products are 
differential magnitudes in the standard Stromgren sys- 
tem. 

During the three consecutive observing seasons be- 
tween 2008 January and 2010 June, the APT acquired 
318 differential brightness measurements of HD 97658 
with respect to the comparison star HD 99518 (V — 7.71, 
B — V = 0.343, F0). We combined the b and y differential 
magnitudes into (b + y)/2 measurements achieving typi- 
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calsingle measurement precision of 1.5-2.0 mmag (|Henrvl 
1x991 . 

The 318 measurements of HD 97658 are plotted in the 
top panel of Figure [6] The second and third observing 
seasons have been normalized to match the mean bright- 
ness of the first season; the second and third year cor- 
rections were 1.75 and 0.70 mmag, respectively. This re- 
moves small year-to-year brightness changes in HD 97658 
and its comparison star and maximizes sensitivity to 
brightness variability on night-to-night timescales. The 
standard deviation of the resulting normalized three-year 
data set is 1.87 mmag, consistent with measurement er- 
ror. Periodogram analysis confirms the absence of peri- 
odic variability between one and 100 days. 

In the second panel of Figure [U the differential mag- 
nitudes are plotted modulo the RV period. Phase 0.0 
corresponds to the predicted time of mid-transit (Ta- 
ble [2]). A least-squares sine fit gives a semi-amplitude of 
0.09 ± 0.14 mmag. This tight limit to photometric vari- 
ability at the RV period supports the hypothesis that 
the period RV signal is due stellar reflex motion from a 
planet in motion. 

The observations near phase 0.0 are replotted on an 
expanded scale in the bottom panel of Figure |U The 
solid curve in the two lower panels approximates the 
depth (0.001 mag) and duration (three hours) of a central 
transit, derived from the orbital elements and assuming a 
water/ice composition for the planet. The uncertainty in 
the time of mid-transit is approximately the width of the 
bottom panel. The vertical error bar in the lower right 
of the transit window corresponds to the ±1.87 mmag 
measurement uncertainty of a single observation. The 
precision and phase coverage of our photometry are in- 
sufficient to detect shallow transits. 

7. DISCUSSION 

We announce two low-mass planets that were reported 
as ano nymous "planet candidates" in Howar d et al.l 
(|2010efl . HD 97658 b is a super-Earth planet with mini- 
mum mass M P sin i = 8. 2±1. 2 M m in a P = 9.494±0.005 d 
orbit around a Kl dwarf star. Gl 785 b is a Neptune-mass 
planet with minimum mass Mp sin i = 21.6 ± 2.0 in 
a P = 74.39 ± 0.12 d orbit also orbiting a Kl dwarf. 

We see no evidence for transits of HD 97658 b, although 
our cphcmcris and photometric phase coverage preclude 
detection of all but the deepest transits of a bloated 
planet. However, given the a priori transit probabil- 
ity of 4%, it is instructive to speculate about the transit 
signatures of vario us possible p l anet c ompositions. Us- 
ing the models in iSeager et al.l (|2007l) , an 8 M$ planet 
composed of pure Fe, MgSiOa, H2O, or H would have 
radii i? p i = 1.3, 1.9, 2.4, and 5.5 Rq, producing transits 
of depth 0.3, 0.6, 1.0, and 5.2 mmag, respectively. These 
homogeneous planet models are oversimplified, but set 
the scale for admixtures of those ingredients. Transits of 
planets made of solids and water would have depths of 
~0. 3-1.0 mmag, while transits of a planet with a signifi- 
cant atmosphere could be much deeper. 

We have no constraints on transits of Gl 785 b be- 
cause the host star is too far South for APT observa- 
tions. The a priori transit probability is 1%. For com- 
paris on, we considered the transiting planets HAT-P- 
11b flBakos et al.ll2010D and HAT-P-26b (lHartman et al.l 
2010), which have masses 26 and 19 M® and radii 
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Fig. 6.— Top panel: The 318 Stromgren (b + y)/2 differential 
magnitudes of HD 97658 plotted against heliocentric Julian Date. 
The standard deviation of these (normalized) observations from 
their mean (dotted line) is 1.99 mmag. Middle panel: The ob- 
servations plotted modulo the RV period. Phase 0.0 corresponds 
to the predicted time of mid-transit. A least-squares sine fit at 
the orbital period yields a semi-amplitude of 0.09 ± 0.14 mmag. 
Bottom panel: The observations near phase 0.0 plotted on an ex- 
panded scale. The duration of a central transit is just three hours 
(±0.0066 phase units); the uncertainty of the transit time is ±0.64 
days (±0.067 phase units). The precision and phase coverage of 
our photometry are insufficient to determine whether or not shal- 
low transits occur. 

4.7 i?^ and 6.3 Rq, respectively. The implied densi- 
ties, 1.38 and 0.42 gem' 3 , suggest that these planets 
have considerable gas fractions. If GL 785 b has a ra- 
dius in the range 4.7-6.3 R®, equatorial transits will be 
4.4-7.8 mmag deep. Such transits would be readily de- 
tectable from the ground, but would require a consider- 
able observational campaign given the transit time un- 
c ertainty of ±2.0 d. 

iFischer fe Valentil (|2005f) showed that the occurrence 
of giant planets with K > 30 ms -1 correlates strongly 
with [Fe/H]. This has been interpreted as support for 
core accretion models of exoplanet formation. However, 
low metallicit y stars might s till be able to form less mas- 
sive planets. IValen ti (2010) noted that stars known to 
host only planets less massive than Neptune (17 Mq) 
tend to be metal poor relative to the Sun. HD 97658 
([Fe/H] = -0.23 ± 0.03, Afsini= 8.2 ± 1.2 M®) and 
G1785 ([Fe/H] = ±0.08 ± 0.03, Msini= 21.6 ± 2.0 M e ) 
are consistent with this tentative threshold. Before in- 
terpreting this physically it is necessary to check for 
metallicity bias in the subsample of stars around which 
sub-Neptune mass planets can be detected with current 
techniques. Further, firmly establishing the apparent 
anti-correlation between host star metallicity and sub- 
Neptune mass planet occurrence is best done with a 
well-controlled sample with uniform detection character- 
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istics, similar to iFischer fc Valentil ([2005) , or with well- 
understood detectability, similar to the Eta-Earth Sur- 
vey. 

We thank the many observers who contributed to 
the measurements reported here. We gratefully ac- 
knowledge the efforts and dedication of the Keck Ob- 
servatory staff, especially Scott Dahm, Hien Tran, and 
Grant Hill for support of HIRES and Greg Wirth for 
support of remote observing. We are grateful to the 
time assignment committees of the University of Cal- 
ifornia, NASA, and NOAO for their generous alloca- 
tions of observing time. Without their long-term com- 
mitment to RV monitoring, these long-period planets 
would likely remain unknown. We acknowledge R. Paul 
Butler and S. S. Vogt for many years of contribut- 
ing to the data presented here. G.W. M. acknowledges 
NASA grant NNX06AH52G. G. W. H. acknowledges sup- 
port from NASA, NSF, Tennessee State University, and 
the State of Tennessee through its Centers of Excellence 
program. This work made use of the SIMBAD database 
(operated at CDS, Strasbourg, France), NASA's As- 
trophysics Data System Bibliographic Services, and the 
NASA Star and Exoplanet Database (NStED). Finally, 
the authors wish to extend special thanks to those of 
Hawai'ian ancestry on whose sacred mountain of Mauna 
Kea we are privileged to be guests. Without their gen- 
erous hospitality, the Keck observations presented herein 
would not have been possible. 

TABLE 4 

Radial Velocities and Shk values for HD 97658 
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Radial Velocity Uncertainty 
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TABLE 5 

Radial Velocities and Shk values for Gl 785 





Radial Velocity 


Uncertainty 




JD - 2440000 


(m s -1 ) 


(ms -1 ) 


Shk 


13237.92941 


1.73 


0.59 


0.2103 


13301.71519 


4.76 


1.13 


0.2260 


13549.02705 


-2.75 


1.02 


0.2040 


13926.01730 


-1.63 


0.56 


0.2023 


13982.83072 


-0.75 


0.50 


0.1963 
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TABLE 5 — Continued 





Rcidicil Velocity 


Uncertainty 




TP) — 9440000 


(m s ) 


l ms ) 
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0.63 


0.1750 
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